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Abstract 

On the basis of a dilatation invariant Lagrangian, governed equations are deter- 
mined for probability density and gauge potential of the non-stationary self-similar 
stochastic system. It is shown that an automodel regime is observed at small time 
interval determined by the Tsallis' parameter q > 1. An exponential falling down 
happens at large time where the dilatation parameter and the partial scale tend to 
constant values. 
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1 Introduction 



Let us consider the one-dimensional random walker with coordinate x at time t - a 
characteristic example of the stochastic system carrying out Levy flights [1]. In the case 
of self-similar system, the corresponding probability distribution P(x, t) is a homogeneous 
function satisfying to the condition 

P{x, t) = a~ a V(K), k = x/a (1) 

where a = a(t) is a time dependent partial scale, k is a dimensionless coordinate, a is a 
self-similarity index. To analyze such a system, it is convenient to use the basic conception 
of Jackson's derivative T> qi whose properties are given in Appendix. Basic advantage of 
this derivative for analyzing self-similar system is that the Jackson's derivative determines 
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rate of the function variation with respect to dilatation q, but not to the shift ckc — > 
as in usual case. This work is devoted to studying a stochastic self-similar system on the 
basis of such a type representation. 

The paper is organized as follows. In Section 2 the governed equations for probability 
density and gauge potential are obtained starting from a dilatation invariant Lagrangian. 
Section 3 deals with the determination of the time dependencies for a characteristic scale 
and probability density. A steady state is shown to realize when the Jackson's derivative 
of the gauge potential equals to zero. This represents a gauge condition, under which 
the probability distribution has the Tsallis' form. With breaking the gauge condition, an 
automodel regime is observed at small time interval determined by the Tsallis' parameter 
q > 1. An exponential falling down happens at large time where the dilatation parameter 
and the partial scale tend to a constant values. Section 4 contains a short conclusions and 
in Section 5 basic properties of the Jackson's derivative are adduced. 



2 Basic equations 

In contrast to the simple case, when the scale a does not depend on the time t, we 
study here a non- stationary self-similar system, for which the value a and the dilatation 
factor q are functions of t. As it is known from the theory of the gauged fields [2], in 
such a case the system is invariant with respect to transformations x — ■> qx, P — > Q q P, 
Qq ~ Q a if the gradient terms V q P, dP/dt are replaced by the elongated derivatives 
(D q + e)P, (d/dt + E)P with dilatational e and temporal E components of the gauge 
potential (hereafter the time t is measured in units of the probability relaxation time). In 
accordance with Eq.(33), it is easy to see that this elongated derivatives are invariant with 
respect to the non-stationary dilatation q — q(t) determined by the follow transformations: 

x — > qx, P — > Q q P] 

(2) 

f^f VqQq (a K V <Q* V * P F^F (1 

e ^ e "^T" (g " 1) ^T^' 

where point denotes the time derivative for brevity. 

Gauge invariant Lagrangian of the corresponding Euclidean field theory is supposed to 
take the form 

C=\[{V q + e)Pf+ l -{Vqef (3) 
where the first term is caused by the gauged dilatation, the second one is the field contri- 
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bution. The respective dissipative function reads 



1 2 



9 

+ 2 




(4) 



where 9 is the relaxation times ratio of the gauge field and probability. As a result, the 
Euler equation 

^ dC dC Of 

v *wd~to = -w z=(p ' e) (5) 



leads to the differential equations with partial derivatives and non-linear terms: 



P + V 2 P=-EP + e 2 P, (6) 
9e + V 2 q e - PV q P = -9EP + eP 2 . (7) 

The first terms in right-hand parts describe dissipative influence of external environment. 
Further, we will take into consideration conserved systems only, so that the time compo- 
nent of the gauge potential being inversely proportional to corresponding relaxation time 
will be put equal to zero (E — 0). 



3 Solution of equations 



In the limits e — > 0, V 2 e — > 0, the obtained equations (6), (7) take the form 



P = -V 2 P, (8) 
9e = PV q P. (9) 

The first of them has the diffusion type but with inverted sign, so that self-similar system 
reveals running away kinetics that is inherent in hierarchical systems [4] . However, such a 
behaviour realizes during short time interval t ~ 9 <ti 1 only. At usual time t > 1, we can 
use the condition 9 1 of adiabatic approximation that will be used everywhere below 
and allow to neglect left-hand side of Eq.(9). As a result, the condition T> q P holds 
true and the system passes to stationary homogeneous regime: 

P(x,t) = const = P st . (10) 
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In much more complicated limit e(t) — > const ^ 0, equation (7) is reduced to the simplest 
form 

V q P = -eP (11) 



and we arrive at the static distribution (10) as before. 

To continue analysis, let us multiply Eq.(6) by factor T> q P and Eq.(7) by T> q t. Then, after 
addition of the obtained results we find 

\{V q Pf + \{V q ef = \{ePf + 1 -C 2 . (12) 

Here we put 

(P - PV q e)V q P = (13) 

and fulfilled integration with constant C 2 /2. As a result, under the assumption 

V q e = -C, C>0 (14) 

and condition V q P ^ 0, the equation (13) arrives at the exponential time dependence 

P oc e~ ct (15) 

whereas the equation (12) is reduced to the form (11). 

The exponential falling-down is known to be not inherent in the self-similar systems [3] 
and, as a consequence, we are need to put C = in Eqs.(12), (14). We arrive then at the 
gauge condition 

V q e = 0, (16) 

according to which the potential e can be a time dependent function but does not vary with 
the system dilatation. Then, equation (6) takes the form P = meaning that the system 
is in steady-state, which probability distribution is obeyed to Eq.(ll). Being accompa- 
nied with Eq.(31), this equation arrives at the condition [a] q = — e. As is ascertained 
in Appendix, the Jackson's q-number [a] q is reduced to the Tsallis' q-logarithm if the 
steady-state probability P st and dilatation q are connected via the follow relation: 



Pr^q a . (17) 
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Then, the above obtained condition gives the Tsallis' distribution [3] 

P st =[l-(q-l)e]^. 



(18) 



With breaking gauge condition (16) the self-similar system gets into non-stationary state, 
whose behaviour is determined by Eq.(6). Accounting definitions (34) arrives this to the 
algebraic form with respect to the Jackson's derivation: 

P = (c 2 -[«]„) P. (19) 

Inserting here Eq.(l) and taking into consideration the condition [a] q = — e and relation 
P = —a~^ 1+a \aV + nV')a (hereafter prime denotes the usual derivation with respect to 
the argument k) we obtain 

a~ x a(aV + kV) = [5a] qq V (20) 

where the factor [5a] qq stands for the term determined by Eqs.(34). 

In the limit q — > oo, one can see with accounting asymptotics (35) that system behaves 
in automodel manner if conditions aq = const, a 3a ~ 4 a = const = t _1 and equation 

kV = (t - a)V (21) 

are implemented. Solution of the equation is V oc K T0 ~ a and the time dependencies of the 
characteristic scale and the probability density read: 

a 3(a-l) = 1 ? p K x r - a t - T> T = T ftt q > : t < r> (22) 

r 3 (a — 1) 

Within the opposite limit g — > 1, a magnitude g in Eq.(20) ought to put time independent 
and we arrive at the long-time dependencies: 

a oc exp (t/r ), P oc x x ° T °~ a exp (-A t), A = ^— |- at g -> 1, t > Aq x . (23) 

The coincidence condition for the time limits r and Ag 1 in dependencies given by Eqs.(22), 
(23) leads to relation 

r = 3(g-l). (24) 
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At last, we consider the case with non-zeroth second Jackson's derivative T> 2 e = [e] qq e 
determined by an index e. Here equation (7) gives the gauge potential 

— ~ [a]q (25) 



1 - [e] qq P~ 2 

that behaves in self-similar manner if the value [e\ qq P~ 2 falls down with g-increase. Sup- 
posing this falling down in power form g~ 7 with positive index 7 — > 0, we obtain the 
needed dependence P(t) oc q T (t), following from Eqs.(22) and condition a(t)q(t) = const, 
if 7 = 2r — 3(e — 1) > 0. As a result, the gauge potential index e is limited by the condition 

e< 1 + - t = l + 2(g-l) (26) 

where the second equality follows from Eq.(24). Under this conditions the equation (6) 
accompanied with approximated result e — [a\ q , following from Eq.(25), arrives at the 
above obtained time dependencies (22). 

The automodel regime is studied to be broken if the factor in right-hand side of Eq.(19) 
e 2 - [a] qq = [a] 2 q [(1 - \s] qq P; 2 )- 2 - l] - [5a] qq = -A (27) 



becomes time independent (here the steady-state probability P st is determined by Eq.(18) 
and we take into account Eqs.(25), (34)). In such a case, the exponential decay (23) is 
characterized by the relaxation time A -1 instead of Aq 1 . Because this regime is inherent 
in small values of q, we can use the limit q — > 1: 



A = A 



a 



e-l 



,2e 



Ao 



a 



(28) 



Finally, under conditions q — > 1, e 2 + (e — l)/(q — 1) «C e 2a when a = [a] q = — e, we have 



A = Ao 



I -2a 



a — 1 



£-1 

9-1 



-2a 



(29) 



4 Discussion 



The above offered formalism is based on the dilatation invariant Lagrangian (3) and dissi- 
pative function (4) to describe the conserved non-stationary self-similar stochastic system. 
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Behaviour of such a system is determined by the probability density distribution (1) and 
the gauge potential e (the latter is reduced to the ratio of the microstate energy to tem- 
perature in usual case of thermodynamic systems). The non-linear differential equations 
with partial derivatives, Eqs.(6), (7) are obtained to analyze the system kinetics. It is 
occurred that under gauge condition (16), when the microstate energy is independent on 
dilatation, the system is in the steady-state characterized by the Tsallis' distribution (18). 
With gauge breaking, the automodel regime (22) realizes for a time less than bounded 
magnitude (24) determined by the difference q — 1. For more values of time, when the 
dilatation becomes constant, the system passes to the usual exponential regime (23). 

Finally, we comment on limitations of our approach. The main of these is that the system 
under consideration is conserved, so that influence of external environment have been 
put equal to zero (the value E = in Eqs.(6), (7)). Accounting such type terms for 
hierarchical systems shows that the time component of the gauge potential is reduced to 
the linear differential operator E = —(d/dx)F(x) + (d 2 /d 2 x)D(x) to be anti-dissipative 
in the physical meaning (here F(x) is a drift force and D(x) is a diffusion coefficient) [5, 
6]. Due to these terms, the above found exponential regime will be suppressed and the 
self-similar system will behave in automodel manner during whole time interval. 



Acknowledgment 

I am grateful to Constantino Tsallis for sending his review [3], which studying inspired 
me to this work. 



5 Appendix. Basic properties of the Jackson's derivative 

The Jackson's derivative is defined by equation 



that is reduced to the usual derivative in the limit q — > 1. Apparently, for a homogeneous 
function with a self-similarity index a the Jackson's derivative is reduced to Jackson's 
q- number [a] q : 



V q f(x) 



/(gap ~ f( x ) 
q-l 



q^l 



(30) 




q a - 1 

q-i 



(31) 
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It is easily to see that the value [a] q — > a in the limit q — > 1 and increases as g"^ 1 
at g — > oo (we propose a > 1). On the other hand, the Tsallis' q-logarithm function 
ln q x = (x q ^ 1 — l)/(q — 1) can be represented in the form of the Jackson's q-number with 
index a — (q — 1) Inx/ lnq. Accompanied Eq.(31) this relation and apparent equality [3] 

ln q (xy) = ln q x + ln q y + (q - 1) (ln q x) (ln g y) (32) 

lead to important rule for the Jackson's derivative: 

V q [f(x)g(x)} = [V q f(x)\ g(x) + f(x) [V q g{x)\ + (q - 1) [V q f(x)] [V q g(x)} . (33) 

The Jackson's derivative of the second order is determined as follows: 



V p V q f(x) = V p {[a] q f(x)} = [a] pq f(x), 

p a — 1 p a \(pq) a — pq\ 

[a} pq = [a] p [a] q + [5a] pq , [a] p = -, [5a] pq = 1 — — — 



(34) 



p-1' L ,Pq ~ (p - l)(pq - 1) 
In proposition a > 1, the value [Sa] pq has the follow asymptotics: 



[Sa] pq -> ^— y at p, g -»• 1, 



(35) 
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